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ABSTRACT 

We present strong detections of methyl cyanide (CH 3 CN), vinyl cyanide 
(CH 2 CHCN), ethyl cyanide (CH 3 CH 2 CN) and cyanodiacetylene (HC4CN) 
molecules with the Green Bank Telescope (GBT) toward the Sgr B2(N) molec- 
ular cloud. Attempts to detect the corresponding isocyanide isomers were only 
successful in the case of methyl isocyanide (CH 3 NC) for its J K = 1 0 - 0 o tran- 
sition, which is the first interstellar report of this line. To determine the spatial 
distribution of CH 3 NC, we used archival Berkeley-Illinois-Maryland Association 
(BIMA) array data for the J K = 4^-3*- ( K = 0 — 3) transitions but no emission 
was detected. From ab initio calculations, the bonding energy difference between 
the cyanide and isocyanide molecules is >8500 cm -1 (>12,000 K). Thus, cyanides 
are the more stable isomers and would likely be formed more preferentially over 
their isocyanide counterparts. That we detect CH 3 NC emission with a single 
antenna (Gaussian beamsize (O# ) = 1 723 arcsec 2 ) but not with an interferometer 
(Qb — 192 arcsec 2 ), strongly suggests that CH 3 NC has a widespread spatial distri- 
bution toward the Sgr B2(N) region. Other investigators have shown that CH 3 CN 
is present both in the LMH hot core of Sgr B2(N) and in the surrounding medium, 
while we have shown that CH 3 NC appears to be deficient in the LMH hot core. 
Thus, large-scale, non-thermal processes in the surrounding medium may account 
for the conversion of CH 3 CN to CH 3 NC while the LMH hot core, which is dom- 
inated by thermal processes, does not produce a significant amount of CH 3 NC. 
Ice analog experiments by other investigators have shown that radiation bom- 
bardment of CH 3 CN can produce CH 3 NC, thus supporting our observations. We 


’NASA Goddard Space Flight Center, Computational and Information Sciences and Technology Office, 
Code 606, Greenbelt, MD 20771 

2 National Research Council Resident Research Associate 

3 Optical Technology Division, National Institute of Standards and Technology, Gaithersburg, MD 20899 
4 National Radio Astronomy Observatory, P. 0. Box 2, Green Bank, WV 24944-0002 


- 2 - 


conclude that isomers separated by such large bonding energy differences are dis- 
tributed in different interstellar environments, making the evaluation of column 
density ratios between such isomers irrelevant unless it can be independently 
shown that these species are co-spatial. 

Subject headings: ISM: abundances - ISM: clouds - ISM: individual (Sagittarius 
B2(N)) - ISM: molecules - radio lines: ISM 


1. INTRODUCTION 

The observations of isomers in the interstellar medium is becoming increasingly more 
common as new and larger molecules are being detected. Currently, 135 molecules have been 
detected in interstellar and circumstellar environments. Excluding diatomic molecules, ~30% 
of all interstellar molecules have observed isomeric counterparts. In addition, as the number 
of atoms in an interstellar molecule increases, so do the number of isomers. For example, 
of the 9 interstellar molecules containing 8 atoms, 5 are isomers. With an empirical for- 
mula C 2 H 4 O 2 , acetic acid (CH 3 COOH), glycolaldehyde (CH 2 OHCHO) and methyl formate 
(CH 3 OCHO) are isomers; and for H 2 C 6 , hexapentaenylidene and triacetylene (HCeH) are 
isomers. Isomerism may provide important clues to the formation of interstellar molecules. 
As such, an important class of isomers are cyanide (-CN) and isocyanide (-NC) species. 

Hydrogen cyanide (HCN), methyl cyanide (CH 3 CN), vinyl cyanide (CH 2 CHCN), and 
ethyl cyanide (CH 3 CH 2 CN) all have large abundances and are easily detectable in a num- 
ber of interstellar environments and are readily seen in the spectra of several molecular line 
surveys (e.g. Friedel et al. 2004, Nummelin et al. 1998). However, while there have been 
detections of HNC in interstellar and cometary environments, there have been no confirmed 
detection of the isocyanides CH 3 NC, CH 2 CHNC, and CH 3 CH 2 NC. The first tentative detec- 
tion of CH 3 NC was toward Sgr B 2 (OH) by Cernicharo et al. (1988) from observations of the 
J = 4 — 3 ,«/ = 5 — 4, and J = 7 — 6 transitions. Though the lines were highly contaminated 
from other known interstellar species, they calculated an CH 3 NC/CH 3 CN abundance ratio 
of ~0. 03-0. 05. Since this initial claim, there have been no further observations to confirm 
the existence of interstellar CH 3 NC, in particular toward the LMH hot core of Sgr B 2 (N), 
the best known source of large molecules. Thus, in this work, we confirmed the existence 
of CH 3 NC and attempted to observe other large isocyanides by observing the same (or sim- 
ilar) low energy transitions of the cyanides CH 3 CN, CH 2 CHCN, CH 3 CH 2 CN, and HC 4 CN 
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and their corresponding isocyanide isomers toward the the high-mass star forming region, 
Sgr B2(N). 


2. OBSERVATIONS AND RESULTS 

The observations using the NRAO 100-m Robert C. Byrd Green Bank Telescope (GBT) 
were made in 2004 February 24, February 29 and March 29. Table 1 lists the molecular line 
parameters of the transitions of the cyanides and isocyanides sought. The species, calculated 
rest frequencies, transition quantum numbers, line strengths, upper level energies, and con- 
tinuum levels are listed in the first six columns of Table 1 . The GBT half-power beamwidth 
(0 B ) is given as 740 "/ v [GHz]. The K-band receiver set used has a frequency range covering 18 
GHz to 22.5 GHz. For the observations reported herein, the average beamwidth < 0 B >=39" 
corresponds to an average Gaussian beamsize of < Sl B >=tv 9 2 b /[A ln(2)] = 1723 arcsec 2 and 
an average beam efficiency < r) B >~0.8. The GBT spectrometer was configured in its 8 
intermediate frequency (IF), 200 MHz, 3-level mode which permits observing four 200 MHz 
frequency bands at a time in two polarizations through the use of offset oscillators in the IF. 
This mode affords 24.4 kHz channel separation. Antenna temperatures are on the T\ scale 
(Ulich & Haas 1976) with estimated 20% uncertainties. The Sgr B2(N-LMH) J2000 pointing 
position employed was a= 17 h 47 m 19 s .8, 8 = -28°22'17" and the LMH systemic LSR source 
velocity of +64 km s _1 was assumed. Automatically updated dynamic pointing and focusing 
corrections were employed based on realtime temperature measurements of the structure 
input to a thermal model of the GBT; zero points were adjusted typically every two hours 
or less using the calibrators 1626-298 and/or 1733-130. Data were taken in the OFF-ON 
position-switching mode, with the OFF position 60' east in azimuth with respect to the ON 
source position. A single observation consisted of two minutes in the OFF source position 
followed by two minutes in the ON source position. The two polarization outputs from the 
spectrometer were averaged in the final data reduction process to improve the signal-to-noise 
ratio. 

Figure 1 shows our attempt to detect cyanide and isocyanide species toward Sgr B2(N) 
with the GBT. The panels in the left column show the detections of CH 3 CN, CH 2 CHCN, 
CH 3 CH 2 CN and HC4CN. We have successfully detected all species in the left column and for 
the first time the Jk u ,k c = 2n — lio transition of interstellar CH3CH2CN is reported. In the 
panels containing the transitions of CH 3 CN, CH 2 CHCN and CH 3 CH 2 CN, the intensity scale 
(ordinate) decreases, showing that the emission from the larger species is in general weaker. 
The abscissa is the radial velocity with respect to the local standard of rest (LSR) calculated 
for the rest frequency of the transition shown at the top left of each panel assuming a source 
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velocity of +64 km s' 1 . Dashed lines show LSR velocities of +64, +73 and +82 km s' 1 . 
Molecules with no hyperfine (HF) structure previously observed with the GBT such as 
propenal (CH 2 CHCHO), propanal (CH 3 CH 2 CHO), and glycolaldehyde (CH 2 OHCHO) show 
absorption components at +64 and +82 km s' 1 and an emission component at +73 km s' 1 . 
In general, the cyanide spectra shown in the left panel of Figure 1 are consistent with 
this source velocity structure. In the cases of CH3CN and CH 2 CHCN, however, the source 
velocity structure is confused by molecular HF splitting. The fiducial on each HF component 
denotes the position corresponding to +64 km s' 1 . Thus, it can be seen that only the F=0-1 
components are clear of any complication. The panels in the right column of Figure 1 show 
the passbands containing the transitions of CH3NC, CH 2 CHNC, CH3CH 2 NC and HC4NC. 
We clearly detect the absorption and emission components at +64, +73 and +82 km s' 1 
from CH3NC and there is possible narrow emission from CH 2 CHNC near +64 km s' 1 . 
The primary absorption component of CH3NC at +64 km s' 1 is more than an order of 
magnitude weaker in intensity than the corresponding feature of CH 3 CN, however this is the 
first interstellar report of this line. No emission or absorption is seen from either CH 3 CH 2 NC 
or HC4NC beyond a 1 a detection limit (~ 2 mK). 


The Jk = 4*-— Zk ( K = 0—3) CH 3 NC observations using the Berkeley-Illinois-Maryland 
Association (BIMA) Millimeter Array 5 were made in 2001 April 20 in its C configuration to- 
ward Sgr B2(N-LMH). The pointing position was the same as used for the GBT. Table 2 lists 
the molecular line parameters of each species in the CH 3 NC passband. The species, the calcu- 
lated rest frequencies, the transition quantum numbers, the product of the line strength and 
the square of the dipole moment (Sfi 2 ), and the upper level energies are listed in the first five 
columns of Table 2. The synthesized beamwidth for these observations is 6 a x0b — 23/ / 9x7/ / l 
corresponding to a Gaussian beamsize of fls=192 arcsec 2 . The quasar 1733-130 was used to 
calibrate the antenna based gains and Uranus was used as the flux density calibrator. The ab- 
solute amplitude calibration of these sources is accurate to within ~20%. The passbands were 
calibrated automatically during data acquisition; a technical description of auto-calibration 
can be found at http://astron.berkeley.edu/~plambeck/technical.html. The spectral win- 
dow containing the CH3NC transitions has a bandwidth of 50 MHz and is divided into 128 
channels for a spectral resolution of 0.39 MHz per channel. All data were combined, imaged, 
and self-calibrated using the MIRIAD software package (Sault, Teuben & Wright 1995). 


Figure 2 shows fiducials for the Jk = 4# — 3*- (K = 0 — 3) transitions of CH3NC in a 
passband dominated by strong emission from CH 3 CH 2 CN which is largely confined to the 
LMH hot core (Liu &: Snyder 1999). The rest frequency located in the upper left of the panel 


5 Operated by the University of California, Berkeley, the University of Illinois, and the University of 
Maryland with support from the National Science Foundation. 
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corresponds to the Jk„,k c = 928 — 827 transition of CH 3 CH 2 CN for an assumed LSR velocity 
of +64 km s -1 . Dashed lines denote LSR velocities of +64, +73 and +82 km s -1 . The 
transition of CH 3 CH 2 CN clearly shows strong emission components at the LMH systemic 
velocity of +64 km s -1 and also at +73 km s" 1 . The Jk — (K = 0 — 3) CH 3 NC 

transitions are labeled in this passband assuming an LSR velocity of +64 km s -i and, except 
for the Jk = 4 3 — 3 3 transition, all lie in a region that is clear from any other molecular line 
emission. However, no transition of CH 3 NC is detected above ala detection limit (~0.3 Jy 
beam -1 ). The 1 a rms noise level fiducial is shown at the left of the panel. 


3. DISCUSSION 

We have presented two disparate sets of observations: GBT single antenna results with 
a Gaussian beamsize of Gb= 1723 arcsec 2 , and BIMA array results obtained with a Gaussian 
beamsize of Cl B =lQ2 arcsec 2 . Given the GBT beamsize is a factor of 9 greater than the 
BIMA beamsize, the spatial scale of the GBT will couple better to the extended emission 
toward Sgr B2(N). Conversely, the BIMA array will tend to couple better to the 28 arcsec 2 
(~5" diameter) of the LMH hot core. Hence, we expect that the GBT observations will 
be characterized by a lower excitation temperature (Tea;) than the excitation temperature 
associated with the BIMA observations. First, we will derive the cyanide and isocyanide 
column densities from the GBT observations. Then, we will show from BIMA observations 
that CH 3 NC must be extended. 

For the GBT observations, following the formalism outlined in Hollis et al. (2004b), the 
total beam-averaged column density of a molecular species obtained with a single element 
radio telescope is given as: 

For emission: 
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In both equations, the line shapes are assumed to be Gaussian; r] B is the telescope 
beam efficiency; T ex is the excitation temperature; AT^AV is the product of the fitted line 
intensity (mK) and line width (km s -1 ); Q r is the rotational partition function; Sfi 2 is the 
product of the transition line strength and the square of the dipole moment (Debye 2 ); and 
E u is the upper rotational energy level (K). In equation (1), v is the transition frequency 
(MHz) and T bg ~2.7 K is the cosmic background temperature. In equation (2), T c is the 
source continuum temperature (K), and Ei is the lower rotational energy level (K). 

Using the GBT at K-band, Hollis et al. (2004b) found a glycolaldehyde excitation tem- 
perature of T ex = 8 K that satisfied the detected integrated line intensities of the absorption 
components at +64 and +82 km s -1 and the emission component at +73 km s -1 . Since the 
LSR velocity components of the cyanide and isocyanide molecules in Figure 1 appear to be 
similar, we used T ex —8 K to find the total column density of each of the detected species. 
Columns (7) to (12) of Table 1 summarize the Gaussian fitting intensities and line widths 
for each detected transition or 3 a upper limits for the intensities of non-detections at LSR 
velocities of +64, +73 and +82 km s -1 . In the cases of CH 3 CN and CH 2 CHCN, multiple 
Gaussian fits were used to estimate the intensities and line widths of the HF components. 
The column densities listed in column (13) of Table 1 are the sum of all detected velocity 
components and, in the cases of CH3CN and CH2CHCN, of all velocity and HF components. 
The column densities of the cyanide molecules are within a factor of 3 or 4 of the CH 3 CN 
column density of (~10 14 cm -2 ). The similarity in column densities between cyanides is also 
seen at higher resolution toward other sources (Remijan et al. 2004b). The column density 
ratio CH3NC/CH3CN is ~0.02 and is relevant if and only if the two species are co-spatial. 
The possible detection of CH 2 CHNC gives a column density ratio CH 2 CHNC/CH 2 CHCN of 
~0.005. 

From Remijan et al. (2003), the total beam averaged column density obtained with an 
interferometer is: 


N t = 2.04 x 
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In equation (3), f Aldv is the product of the fitted line intensity (Jy beam -1 ) and line width 
(km s -1 ), 6 a and 9 b are the FWHM beam widths (arcsec) and B = + D s ] is the 

beam filling factor where fls is the solid angle subtended by the synthesized beam of the 
interferometer and is the solid angle subtended by the source emission. For the BIMA 
array data, 9 a x9 b = 23."9x7."l. Columns (7) to (12) of Table 2 summarize the Gaussian 
fitting intensities and line widths for CH 3 CH 2 CN and the 3 a upper limits for the intensities 
of the non-detections at LSR velocities of +64, +73 and +82 km s -1 . The column density 
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of CH3CH2CN in column (13) is based on the excitation temperature of methanol (CH3OH) 
of 170 K (Pei, Liu & Snyder 2000). 


No emission for the J = 4 — 3 (K = 0 — 3) transition of CH 3 NC was seen in the 
BIMA array SgrB2(N) data. However, if CH3NC and CH 3 CN are co-spatial (i.e., our GBT 
CH3NC/CH3CN ratio of ~0.02 is valid), then observations of CH 3 CN by other investigators 
which yield large CH3CN column densities predict that we should have detected CH 3 NC 
with the BIMA Array if the source is compact. In an extensive survey encompassing the 
J = 12 — ♦ 11, J = 13 — ► 12 and J = 14 — > 13 transitions of CH 3 CN at 1 mm wavelengths, 
Nummelin et al. (2000) obtained a CH 3 CN beam averaged column density of ~4.2xl0 15 cm -2 
for a 352 K LMH hot core temperature. These observations were performed using the 15m 
SEST telescope with a ~600 arcsec 2 beamsize. In addition to the the column density and 
temperature, Nummelin et al. also found a beam filling factor of the CH 3 CN emission toward 
the Sgr B2(N) hot core of 0.014. In their data analysis, the column density, rotational 
temperature and beam filling factor were given as free parameters. From these values, an 
integrated line intensity was calculated for each transition and was compared to the observed 
integrated line intensity until a “best-fit” model was obtained. The best-fit model was found 
by finding the minimum of the reduced x 2 function (Nummelin et al. 2000). Assuming the 
model explains the physics accurately, the measured values of the rotational temperature, 
column density and beam filling factor by Nummelin at al. toward the Sgr B2(N) hot core are 
reported at a 68.3% confidence interval (1 a). Given the ~600 arcsec 2 beamsize of the 15m 
SEST telescope, and a beam filling factor of 0.014, the measured source size of the CH3CN 
emitting region toward the Sgr B2(N) hot core is ~8.5 arcsec 2 . This size is very similar 
to the CH3CN cores toward W51el (7.9(6)arcsec 2 ) and W51e2 (6.0(9)arcsec 2 ) (Remijan 
et al. 2004a). Furthermore, the size is also similar to the emitting regions of CH 2 CHCN 
and CH 3 CH 2 CN towards the Sgr B2(N) hot core measured by Liu & Snyder (1999) at 
subarcsecond resolution. Thus, ~8.5 arcsec 2 is a reasonable measure of the CH 3 CN emitting 
region toward the Sgr B2(N) hot core. Using the measured rotational temperature of 352 K, 
the scaled Nummelin et al. (2000) CH 3 CN column density within the 192 arcsec 2 beamsize 
of the BIMA array is ~9.9x 10 16 cm -2 . Using the column density ratio of 0.02, the expected 
CH 3 NC column density should be ~2xl0 15 cm -2 for co-spatial isomers. 


Solving equation (3) for J A Idv assuming a column density of.~2xl0 15 cm -2 , we should 
have detected a CH3NC integrated line intensity of ~57 Jy beam -1 km s -1 for the Jk = 
4 0 — 3 0 transition. Assuming a typical line width of 15 km s -1 for the LMH hot core 
of Sgr B2(N), the expected peak intensity of the 4 0 — 3o transition is ~3.8 Jy beam -1 . 
However, the BIMA array data shows no CH 3 NC emission greater than the 1 a rms upper 
limit of ~0.3 Jy beam -1 . The most likely explanation for the lack of a CH 3 NC detection 
is that the emission is extended and the array cannot couple to it. In the case of extended 
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emission, the ratio of the measured line intensity, or in our case the 1 a rms upper limit of 
~0.3 Jy beam -1 , to the expected line intensity of ~3.8 Jy beam -1 is equal to the ratio of 
the solid angle subtended by the synthesized beam of the interferometer (Q B ) to the solid 
angle subtended by the source emission (f }$)■ This is because the synthesized beam of the 
interferometer is only collecting a fraction of the flux from the extended emission region which 
we assume to be uniformly distributed. Therefore our observations indicate that < 

0.3/3. 8. Solving for the source size assuming a Gaussian beamsize of 192 arcsec 2 , we find 
fls>2432 arcsec 2 or equivalently, a linear size of >100 /, (1.7'), which is much larger than the 
typical LMH hot core diameter of ~5" (Hollis et al. 2003). This strongly suggests that the 
spatial distribution of CH3NC is extended and would not show any significant concentration 
toward the LMH hot core. 

DeFrees, McLean, & Herbst (1985) theoretically calculated a CH3NC/CH3CN abun- 
dance ratio for dense interstellar clouds in the range of 0. 1-0.4 by assuming that both species 
are formed from ion precursors and destroyed by gas phase, ion molecule reactions. The 
DeFrees et al. calculations show that when CH3 + collides with HCN, a collision complex is 
formed that rapidly equilibrates to protonated methyl isocyanide (CH 3 NCH + ) and proto- 
nated methyl cyanide (CH 3 CNH + ) due to the large amount of energy produced from the 
initial collision. As the ensemble of ions continues to relax, the internal energy of CH3CNH + 
becomes insufficient to overcome the barrier to isomerize to CK 3 NCII + . The end result finds 
CH3CN formed preferentially over CH3NC. We performed ab initio calculations in order to 
quantify the energy differences between cyanide and isocyanide isomers in this work. Table 
3 lists each molecular species (col. 1), the zero-point corrected bonding energy (ZPCBE in 
col. 2) and the relative ZPCBE (col. 3). Column (2) ZPCBE represents the complete disso- 
ciation of the molecule into its constituent atoms. Note that the relative bonding energies 
between isomers axe very large (>8500 cm -1 ). Thus, the cyanide molecules are much more 
stable than their isocyanide isomers, and therefore it is highly unlikely that any thermal 
process can satisfactorily explain the formation of isocyanide molecules. This indicates that 
a non-thermal process is necessary to convert a cyanide into an isocyanide species. 

In recent laboratory experiments performed by Hudson & Moore (2004), cyanide species 
were bombarded by both protons and UV-photons to determine what would happen to the 
structure of the molecule. Bombarding a pure ice sample of CH 3 CN at 24 K, the proton 
experiment yielded ketenimine (CH 2 CNH), another isomer of CH3CN, as well as CH 3 NC in 
a slightly lower abundance. Similarly, the UV-photon experiment at 12 K yielded CH 3 NC 
as the primary product. These products were also seen when CH 3 CN was embedded in N 2 
ice and then irradiated. Hudson & Moore (2004) also note that the temperature in both 
experiments did not affect the formation of the end products. In an UV-photon experiment 
performed at 100 K, the destruction of CH3CN was faster than at 12 K, but the same end 
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products were seen. This suggests that the reactions were driven more by the radiation 
dose than the sample temperature (Hudson & Moore 2004), reinforcing our suggestion that 
non-thermal processes may be the primary route to the formation of interstellar isocyanides. 
Laboratory experiments involving CH 2 CHCN, CH 3 CH 2 CN, and cyanoacetylene (HC 2 CN), 
all showed similar results for isocyanide isomers. It is likely that a similar reaction involving 
HC4CN will produce HC4NC. Furthermore, these laboratory experiments also suggest that 
all isocyanides are likely to be formed under the conditions which interstellar CH 3 NC was 
detected. Based on the column density ratios and upper limits in Table 1, we did not de- 
tect any emission from CH 3 CH 2 NC and HC 4 NC because of the possibility that we did not 
integrate long enough to detect these species. On the other hand, Hudson &: Moore (2004) 
also conducted experiments with cyanides embedded in H 2 0 ice and found that hydrogen 
addition and subtraction reactions were taking place in both radiation bombardment reac- 
tions in which no isocyanides were formed. In these experiments HC 2 CN hydrogenated to 
CH 2 CHCN and then to CH 3 CH 2 CN. Based on the aforementioned laboratory results, we 
will suggest how large cyanides and their corresponding isocyanide isomers may be produced 
in the extended regions surrounding Sgr B2(N). 


The Sgr B2(N) region is known to contain widespread shocks (Chengular & Kanekar 
2003). Such shocks have been presumed to be responsible for the formation and distribution 
of large aldehyde molecules toward this region (Hollis et al. 2004a, b). Furthermore, it is 
also well known that there is ambient UV flux and cosmic ray flux in these regions. On 
the surfaces and mantles of dust grains, larger species such as CH 3 CN and HC 2 CN may 
be present. CH 2 CHCN and CH 3 CH 2 CN are formed as a result of H-addition reactions, 
thus depleting HC 2 CN. As a shock passes, these large molecules are expelled from the grain 
surface into the gas phase. At this point UV and cosmic rays may interact with the cyanide 
species and form isocyanides. As a result, the isocyanides show an extended distribution. In 
hot cores, there is a lower UV and cosmic ray flux due to the large column density of gas 
and dust. Furthermore, thermal processes are the primary mechanism to liberate molecules 
off grain surfaces (Tielens & Hagen 1982; Hasegawa, Herbst, & Lueng 1992). Indeed, this 
is initially why these regions were believed to be the ideal places to form large molecules 
(Ehrenfreund & Charnley 2000). Also, in these compact regions, it is believed warm gas 
phase chemistry in addition to grain surface chemistry can form large molecules which may 
account for the enhanced abundances of CH 2 CHCN and CH 3 CH 2 CN in hot core regions (Liu 
& Snyder 1999, see also CH 3 CH 2 CN column density in Table 2 with respect to Table 1). 
However, thermal processes are not enough to form isocyanide species after the cyanides have 
been expelled from the grain and hence there should be a depletion of isocyanides towards 
hot cores. Our observations support this hypothesis. 


In summary, we have calculated the chemical bonding energies for cyanide molecules 
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studied herein (CH 3 CN, CH 2 CHCN, CH 3 CH 2 CN, & HC 4 CN) and found they are much 
greater (>8500 cm -1 ) than that of the corresponding isocyanide isomers. As a consequence, 
cyanide interstellar molecules are most likely preferentially formed at the expense of iso- 
cyanide isomers. Moreover, in the case of CH 3 CN, it is seen in the LMH hot core and in 
the surrounding medium while CH 3 NC appears to be deficient in the LMH hot core. This 
suggests that large-scale, non-thermal processes such as shocks or enhanced UV flux in the 
surrounding medium account for the conversion of CH 3 CN to CH 3 NC while the LMH hot core 
is dominated by thermal processes which cannot produce a significant amount of CH 3 NC. 
Hudson & Moore (2004) have shown that radiation bombardment of CH 3 CN in ice analog 
experiments can produce CH 3 NC, thus tending to support what is observed. Additionally, 
since CH 3 NC is deficient in the LMH and CH 3 CN is prevalent, this suggests that the relative 
abundance ratio between these two isomers is not meaningful because the molecules are not 
co-spatial. This may well apply to other cyanide and isocyanide isomers. 

J.M.H. gratefully acknowledges research support from H.A. Thronson, Assistant Asso- 
ciate Administrator for Technology, NASA Science Mission Directorate. 
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Table 1. GBT Observations - Cyanide and Isocyanide Molecular Line Parameters and Column 

Densities 


Species 

Frequency 

Transition 

S 

E„ 

T c 

(+64 km 

8" 1 ) 

(+73 km 

s’ 1 ) 

(+82 km 

s’ 1 ) 

Nx 

N t (XNC) 







<*TJ) 

(AV) 

(AH) 

(AV) 

(AT-) 

(AV) 

X10~ 13 

N t (XCN) 


(MHz) 



(K) 

(K) 

(mK) 

(km/s) 

(mK) 

(km/s) 

(mK) 

(km/s) 

(cm -2 ) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

CH 3 CN a 

18396.7252(7) 

lo-0o 
F = 1-1 

1.332 

0.88 

10.6 

-275(25) 

5.8(7) 

160(20) 

7.9(9) 

<7 





18397.9965(6) 

lo-Oo 
F = 2 - 1 

2.220 

0.88 

10.6 

-400(15) 

7.0(3) 

310(19) 

5.9(4) 

<7 


10.8(20) 



18399.8924(3) 

lo-Oo 
F = 0- 1 

0.444 

0.88 

10.6 

-180(25) 

8.3(9) 

<7 


<7 



0.02(2) 

CH 3 NC b 

20105.754(1) 

lo-0o 

4.000 

0.88 

10.2 

-48(2) 

8.4(4) 

11(3) 

6.4(9) 

<7 


0.2(1) 


CH 2 CHCN c 

18512.158(13) 

2i2-lll 

F = 2- 1 

0.375 

3.51 

11.1 

-62(2) 

4.9(1) 

<7 


<7 





18513.311(7) 

2i2-ln 
F = 3 - 2 

0.701 

3.51 

11.1 

-82(2) 

4.9(1) 

<7 


<7 


37.2(30) 



18514.400(27) 

2i2-ln 
F = 1 - 0 

0.167 

3.51 

11.1 

-19(3) 

4.1(5) 

<7 


<7 



0.005(2) 

CH 2 CHNC d 

19993.916(1) 

2i2-lii 

1.500 

3.67 

10.5 

8(2) 

1.5(4) 

< V 


<7 


n o/k\ 
v'.~v v / 


ch 3 ch 2 cn® 

18377.721(3) 

2h-1io 

1.500 

1.56 

10.0 

-24(2) 

2.8(2) 

20(2) 

2.8(2) 

-20(2) 

2.8(2) 

3.5(4) 

<0.3 

CH 3 CH 2 NC f 

18803.209(10) 

2i2-ln 

1.500 

1.56 

9.8 

<7 


<7 


<7 


<1.1 


HC 4 CN« 

18638.617(2) 

7-6 

7.000 

2.70 

10.9 

-63(2) 

5.8(2) 

97(1) 

16.0(5) 

-59(2) 

8.0(3) 

6.1(2) 

<0.03 

HC 4 NC h 

19616.504(2) 

7-6 

7.000 

2.84 

9.7 

<5 


<5 


<5 


<0.2 



Note. — Uncertainties in the frequencies refer to the least significant digit (Taylor Kuyatt 1994) and are 2 er values (coverage factor 2). Upper limits of the 
line intensities are 3 <x values. All other listed uncertainties are 1 a. 

“Rest Frequencies from Boucher et al. 1980; =3.922(1) from Gadhi et al. 1984 

b Rest Frequencies from Kukolich 1972; /i tt =3. 830(60) from Ghosh et al. 1953 

c Rest Frequencies from Gerry et al. 1979; n a =3.815(12), /*& =0.894(68) from Stolze, M. & Sutter D. H. 1985 

d Rest Frequencies from Yamada & Winnewisser 1975; #x a =3.470(30), ^{,=0.790(100) from Bolton, Owen, &: Sheridan 1970 

e Rest Frequencies, /n a =3.850(10) and = 1.230(60) from Lovas 1982 

f Rest Frequencies from Fliege & Driezler 1985; n a =3. 790(20) and Pb=l. 310(40) from Anderson & Gwinn 1976 
g Rest Frequencies from Winnewisser et al. 1982; /i 0 =4.330(30) from Alexander et al. 1976 
h Rest Frequencies from Botschwina et al. 1998; =3.590 calculated in this work 
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Table 2. BIMA Array Observations - Molecular Line Parameters and Column Densities 


Species 

Frequency 

Transition 

S M 2 

E„ 

(+64 km 

® 1 ) 

(+73 km 

8 *) 

(+82 km s *) 

N t 






rs Tj 

(AV) 

(AT) 

(AV) 

(A/) 

<AV) 



(MHz) 


(Debye 2 ) 

(K) 

(J/bm) 

(km/s) 

(J/bm) 

(km/ 8 ) 

(J/bm) 

(km/s) 

(cm -2 ) 

(1) 

( 2 ) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

( 10 ) 

(ID 

( 12 ) 

H 2 CCCC a 

80383.96(2) 

9oo-8o8 

153.8 

19.3 

<0.3 


<0.3 


<0.3 


< 6 . 0 xl 0 14 

HCOOCH 3 b 

80395.135(126) 

928 -909 E 

0.7 

28.8 

<0.3 


<0.3 


<0.3 


<4.7X 10 17 

CH 3 CH 2 CN 

80404.898(12) 

928-8 2 7 

126.1 

23.8 

4.6(1) 

16.0(8) 

1.7(3) 

5.1(6) 

0 . 6 ( 2 ) 

13.6(4.0) 

3.1X10 16 

CH 3 NC 

80405.504(120) 

4 3 -3 3 

29.6 

73.0 

<0.3 


<0.3 


<0.3 



CH 3 NC 

80414.660(120) 

4 2 -3 2 

50.7 

37.4 

<0.3 


<0.3 


<0.3 



CH 3 NC 

80420.062(120) 

4i-3i 

63.3 

16.0 

<0.3 


<0.3 


<0.3 



CH 3 NC c 

80421.910(120) 

4q‘3 0 

67.6 

8.9 

<0.3 


<0.3 


<0.3 


<4. Ox 10 13 


*Rest Frequency and Line Strength (S/i 2 ) from Killian et al. 1990 
b Rest Frequency and Line Strength (S** 2 ) from Bauder 1979 
c N'r calculated using the strongest expected component. 


Table 3. Chemical Bonding 
Energy 


Species 

(i) 

ZPCBE a 

(cm -1 ) 

( 2 ) 

Relative 

ZPCBE 

(cm -1 ) 

(3) 

CH 3 CN 

-29,070,027 

0 

CH 3 NC 

-29,060,541 

9,486 

CH 2 CHCN 

-37,408,030 

0 

CH 2 CHNC 

-37,399,372 

8,658 

CH 3 CH 2 CN 

-37,672,368 

0 

CH 3 CH 2 NC 

-37,663,671 

8,697 

HC 4 CN 

-53,818,766 

0 

HC 4 NC 

-53,807,547 

11,219 


Note. — The Chemical Bonding Energies were 
Obtained at the MP2/aug-cc-pVTZ Level of The- 
ory (Frisch et al. 2004; Dunning 1989; Kendall, 
Dunning & Harrison 1992; Woon & Dunning, 
1993) 

a Zero- Point Corrected Bonding Energy 
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Fig. 1. — GBT cyanide and isocyanide spectra toward SgrB2(N) at 24.4 kHz channel spacing. 
Transition quantum numbers are shown in each panel. Each spectrum was processed with 
a median filter to remove instrumental slopes in the bandpass. The abscissa is the radial 
velocity with respect to the LSR calculated for the rest frequency of the transition shown 
(see Table 1) at an assumed source velocity of +64 km s -1 . Dashed lines show LSR velocities 
at +64, +73 and +82 km s -1 . In the case of CH 3 CN and CH 2 CHCN the HF structure is 
also seen and labeled for an LSR velocity of +64 km s -1 . 

Fig. 2. — BIMA array spectral passband toward the LMH hot core of Sgr B2(N) containing 
the Jk — 4# — 3k (K = 0 — 3) transitions of CH 3 NC. In this case, the rest frequency 
corresponds to the CH 3 CH 2 CN transition at 80.4049 GHz for a LSR velocity of +64 km s -1 . 
However, there are also emission components at +73 and +82 km s -1 . Finally, the locations 
of the transitions of CH 3 NC, HCOOCH 3 and H 2 CCCC (Table 2) are also labeled assuming 
an LSR velocity of +64 km s -1 . 
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Figure 2. 





